Octonionic electromagnetic and gravitational interactions and dark matter 



Zihua Wen£] 

School of Physics and Mechanical & Electrical Engineering, 
Xiamen University, Xiamen 361005, China 

Based on the Maxwellian quaternionic electromagnetic theory, the electromagnetic interaction, 
the gravitational interaction and their coupling influence with the dark matter field in the octonionic 
space are discussed. The paper reveals the close relationships of dark matter field with electromag- 
netic field and gravitational field. In the dark matter field, the research discloses the movement 
properties of dark matter in certain conditions and its influence on ordinary matter movements. In 
the electromagnetic field, the variation of field energy density has direct effect on the force of electric 
charge and current. In the modified gravitational field, the near coplanarity, near circularity and 
corevolving of planetary orbits are deduced from the equations set. The change of centrifugal force 
of celestial body leads to fluctuation of revolution speed, when the field energy density, celestial 
body's angular momentum in its sense of revolution or its mass spatial distribution varied. The 
results explain that some observed abnormal phenomena about celestial bodies are caused by either 
modified gravitational field or dark matter. 
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I. INTRODUCTION 

Nowadays, there still exist some movement phenomena of celestial body which can't be explained by current 
gravitational theory. Therefore, some scientists doubt the universality of current gravitational theory, and then bring 
forward the modified gravitational theory and hypothesis of dark matter to explain the phenomena of abnormal 
rotation velocity of galaxy and the associated 'missing mass' etc. The hypothesis of dark matter believes that, there 
exist not only various gleamy ordinary matters, but also one kind of mysterious matter which exerts the gravitational 
pull, but neither emits nor absorbs light [![. This kind of matter is known as dark matter, which owns mass and has 
effect on ordinary matter movements 0. Most scientists believe in its existence and multiformity in the universe. 

A new insight on the problem of the dark matter can be given by the concept of the octonion space. The under- 
standing of time and space experiences an evolutional period 0]. In 17th century, I. Newton thought the time and 
space in nature were separate. And the time and space were both independent of matter. In 20th century, A. Einstein 
believed that, the time and space in nature were integrated and commonly known as the spacetime. The spacetimc 
related to the movement status and matter. The precursors include R. Descartes and A. Einstein etc. held the idea 
that space and time were the extension of matter. And there was no such spacetime without matter and hence no 
empty spacetime Q. 

In the paper, we believe that electromagnetic and gravitational interactions are different, equal and relative. And 
they possess different and four-dimensional spacetimes respectively. These two types of different and four-dimensional 
spacetimes are equal and independent, and can not be superposed directly. The viewpoint can be summarized as 
'SpaceTime Equality Postulation', that is: each fundamental interaction possesses its own spacetime, and all these 
spacetimes are equal. 

According to previous research results and the 'SpaceTime Equality Postulation', the electromagnetic and gravita- 
tional interactions can be described by their four-dimensional spacetimes. Based on the conception of space verticality 
etc., two types of four-dimensional spacetimes can be united into an eight-dimensional spacetime. In eight-dimensional 
spacetime, the electromagnetic interaction and gravitational interaction can be equally described. So the character- 
istics of electromagnetic and gravitational interactions can be described by a single eight-dimensional spacetime 
uniformly. 

The paper modifies the gravitational theory of ordinary matter and dark matter, and draws some conclusions 
which are consistent with the Maxwellian electromagnetic theory and the Newtonian gravitational theory. A few 
prediction associated with movement feature of dark matter can be deduced, and some new particles can be used for 
the candidate of dark matter. 
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II. ELECTROMAGNETIC FIELD AND QUATERNION SPACE 

The electromagnetic theory can be described with the quaternion algebra In the treatise on electromagnetic 
theory, the quaternion algebra was first used by J. C. Maxwell to describe the equations set and various properties of 
the electromagnetic field @. The spacetime, which is associated with electromagnetic interaction and possesses the 
physics content, is adopted by the quaternion space. 

In the quaternion space of electromagnetic interaction, the radius vector K = (ro, r 1; r 2 , r 3 ) , and the base E = (1, 
*ii *2, i 3 )- Where, ro = ct, c is the speed of light beam, and t denotes the time. 

The quaternion differential operator and its conjugate operator Q* are defined as 

= d + iidi + fa.d 2 + i 3 d 3 , 0* = d - fadi - i 2 d 2 - i 3 d 3 . (1) 

where, di — d/dri, i = 0, 1, 2, 3 . The mark (*) represents the quaternion conjugate. 
The electromagnetic field potential A(ao, flti, a 2l 0,3) is defined as 

A = 0°X = a + i\a\ + h a 2 + ^303 (2) 
where, X is the physical quantity. 



A. Electromagnetic field strength 

The electromagnetic field strength B(&o, &i, b 2 , b 3 ) is defined as 

B = o A 

= b 00 + (*1&01 + hb 2 + «3^03) + (*1&32 + «2&13 + ^^l) 

= (d cio - diai - d 2 a 2 - d 3 a 3 ) 

+ {ii(d ai + dia ) + i 2 (d a 2 + d 2 a ) + i 3 {d Q a 3 + d 3 a )} 

+ {k{d 2 a 3 - d 3 a 2 ) + i 2 (d 3 ax - 9 x a 3 ) + i 3 {dia 2 - d 2 ai)} (3) 

where, the first term in the right of equal mark is gauge definition, the second term is electric field intensity (601 , ^02 , ^03) 
and the third term is flux density {b 32 , b\ 3l b 2 \). Taking 60 = boo = 0, then the gauge condition of field potential in 
electromagnetic field can be achieved. 



TABLE I: The quaternion multiplication table. 

1 w h i3_ 

1 1 ii 12 fa 

ii ii —1 fa —fa 

fa fa —fa —1 fa 

fa fa fa —fa — 1 



B. Electromagnetic field source 

The electromagnetic field source S(so, si, s 2 , s 3 ) is defined as 

-fiS= 0* o B 

= (d b + dih + d 2 b 2 + d 3 b 3 ) 

+ {iiidoh - dibo) + i 2 (d b 2 - d 2 b ) + i 3 {d b 3 - d 3 b )} 

+ {fa(d 3 b 2 - d 2 b 3 ) + i 2 {dib 3 - d 3 h) + i 3 {d 2 bi - <9i& 2 )} 
= 2 (a + faai + i 2 a 2 + i 3 a 3 ) (4) 

where, 2 = 0* • H is a coefficient. 60 = . Taking § = (so + iiSi + i2S 2 + i 3 s 3 ), then equations set of field 
potential of the electromagnetic field can be attained 

2 (a + faai + i 2 a 2 + i 3 a 3 ) = -fi(s + fasi + i 2 s 2 + i 3 s 3 ) (5) 
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TABLE II: Equivalent representation of the electromagnetic field equation. 





traditional representation 


equivalent representation 


field potential 


A 








¥> 


field strength 


B = Vxi 


B = Vxi 




E = -(v</3 + dA/dt) 


E' = S7<p + dA/dt 




D — e E 


D — e E 




H = B 1 fi 


H = 


gauge condition 


7- A + d(ip/c)/d(ct) = 


X7-A - d(ip/c)/d(ct) = 


field source 


V 2 w - d 2 io/d(ct) 2 - -o/e 


V 2 m + d 2 u>/d(ct) 2 — —ole 




si 2 A - d 2 A/d(ct) 2 = -nj 


v 2 ^ + a 2 A/a(ci) 2 = -/xj 


Maxwell equation 


V x H=dD/dt + j 


V x H= dD'/dt + j 




V x E = -dB/dt 


V x E" = dB/dt 




V • D= p 


V • D'=-p 




V • B = 





C. Variation rate of field source 



The variation rate NJoijiiJZija) of electromagnetic field source is defined as 

J= 0*o§ 

= (doso + disi + d 2 s 2 + d 3 s 3 ) 

+ {ii(dasi - diso) + i2(d s 2 - d 2 s ) + i 3 (dos 3 - d 3 s )} 

+ {ii{d 3 s 2 - d 2 s 3 ) + i 2 (d!S 3 - d 3 s x ) + i 3 (d 2 s 1 ~ <9is 2 )} (6) 

Taking the gauge condition jo = 0, then the continuity equation of field source of the electromagnetic field can be 
achieved. 

In the quaternion space, various properties among the field potential, field strength and field source of the elec- 
tromagnetic interaction can be described by the quaternion. Those conclusion and results are similar to that of 
equivalent representation of the electromagnetic field in Table 2. According to the 'SpaceTime Equality Postulation', 
it can be deduced that the spacetime derived from gravitational interaction is supposed to be the quaternion space 
also. Given that the coupling influence of fundamental interactions is neglected, the gravitational interaction is similar 
to electromagnetic interaction and can also be described by quaternion. 



III. ELECTROMAGNETIC-GRAVITATIONAL FIELD AND OCTONION SPACE 



The electromagnetic and gravitational interactions are interconnected, unified and equal. Both of them can be 
described in the quaternion space. By means of the conception of space expansion etc., two types of quaternion 
spaces can combine into an octonion space [7|. In the octonion space, various characteristics of electromagnetic and 
gravitational interactions can be described uniformly, and some correlative equations set can be obtained. 



quaternion 




equivalent 




traditional 


electromagnetic 




electromagnetic 




electromagnetic 


field theory 


< ► 


field theory 


< ► 


field theory 



FIG. 1: Equivalent description in the electromagnetic field theory. In the quaternion space, the various properties 
among field potential, field strength and field source of electromagnetic interaction can be described by quaternion. 
Those conclusion and results are similar to that of equivalent representation of the electromagnetic field. 
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The base E 9 = (1, i\, i^, i 3 ) of quaternion space of gravitational interaction (G space, for short) is independent of 
the base E e of quaternion space of electromagnetic interaction (E space, for short). Selecting E e = (1, i 1; i 2 , i 3 ) o I 
= (Io, Ii, I2, I3). And then the bases E g and E e constitute a single base E of the octonion space. 

E = E g +E e = {l,i u V2,i3,Io,Iul2,l3) (7) 

The radius vector R(ro, ri, r 2 , r 3 , Rq, R 2 , R3) in the octonion space is 

K = (r + kn + hr 2 + i 3 r 3 ) + (IqRq + hRi + I2R2 + I3R3) (8) 

where, tq = vot, i?o = VoT. vq = Vo = c is the speed of light beam, t and T denote the time. 
The octonion algebra is the alternative algebra, so the octonions Qi and Q 2 satisfy 

Qi o (Qi o Q 2 ) = (Qi o Q x ) o Q 2 , Q x o (Q 2 o Q 2 ) = (Q 1 o Q 2 ) o Q 2 . 

The octonion differential operator and its conjugate operator 0* are defined as 



= 9 + e , 0* = 0; + 0: ; (9a) 

9 = d g0 + i\dgi + i 2 d g2 + i 3 d g3 , 0* = d g0 - iid gl - i 2 9 g2 - i 3 d g3 ; (9b) 

e = Iodeo + hdei + hd e 2 + hd e 3 , 0* = - Iod e o - hdei - hde2 - hd e3 . (9c) 

where, d gi — d/dri ; d ei — d/dRi ; i = 0, 1, 2, 3 . The mark (*) represents the octonion conjugate. 

The field potential A(oo, ai,a 2 ,a 3 ,A , Ai,A 2 ,A 3 ) in the electromagnetic-gravitational field is defined as 

A = o X = (a + iiai + i 2 a 2 + i 3 a 3 ) + k a (I A + I X A X + I 2 A 2 + I 3 A 3 ) (10) 



where, X = X s + fc x X e ; X fl and X e are the physical quantity in G space and E space respectively; A = A g + k a A e : 
A g = (do, oi, a 2 , a 3 ) and A e = (A , A\, A 2l A 3 ) are the field potential in G space and E space respectively; k a and k x 
are coefficients. 



A. Field strength 

The field strength B(6 , bi,b 2 , b 3 , B 0} Bi,B 2 ,B 3 ) is defined as 

B = o A 

= (6 + kbi + i 2 b 2 + i 3 b 3 ) + k b (I B a + I X B X + hB 2 + I 3 B 3 ) 
= [{d g oa a - dgiai - d g2 a 2 - d g3 a 3 ) 

+ii(d gQ a 1 + d g ia ) + ia{d g oa 2 + d g2 a ) + i 3 (d g0 a 3 + d g3 a ) 

+i\(d g2 a 3 - d g3 a 2 ) + i 2 {d g3 ai - d gl a 3 ) + i 3 (d gl a 2 - d g2 ai)] 

+ [I (d e0 a + d e iai + d e2 a 2 + d e3 a 3 ) 

+Ii(d el a - d e0 ai) + I 2 (d e2 a - d e0 a 2 ) + I 3 (d e3 a - d eQ a 3 ) 
+Ii(d e3 a 2 - d e2 a 3 ) + I 2 {d el a 3 - d e3 ai) + 7 3 (<9 e2 ai - d e ia 2 )] 
+k a [I (dg Ao - d g iA-i - d g2 A 2 - d g3 A 3 ) 

+7i(a 9 iA + dgoAi) + I 2 (d g2 A + d ga A 2 ) + I 3 (d g3 A + d g0 A 3 ) 
+h{d g3 A 2 - d g2 A 3 ) + I 2 {d gl A 3 - dg^) + htdgiAt - d gl A 2 )] 
+k a [-(d ea A + d eX A x + d e2 A 2 + d e3 A 3 ) 
+i 1 (d e0 A 1 - d el A ) + i 2 (deaA 2 - d e2 A Q ) + i 3 (d e0 A 3 - d e3 A ) 

+ il(d e3 A 2 - e2 A 3 ) + h(delA 3 - de^) + Z 3 (9e2^1 - d el A 2 )} (11) 

where, the first two terms are the field strength of the gravitational interaction, the last two terms are the field 
strength of the electromagnetic interaction. B = R g + k{,M e ; M g = (bi,b 2 ,b 3 ) and B e = (B\, B 2 , B 3 ) are the field 
strength in G space and E space respectively; kb is a coefficient. Selecting gauge equations b = and B — can 
simplify definition of field strength in the octonion space. 
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TABLE III: The octonion multiplication table. 
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Jo 


-is 


-ii 
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-1 



B. Field source 

The field source §(so, si, s 2 , S3, So, Si, S 2 , S 3 ) is defined as 

- fiS = 0* o B 

= [{dgob + d gl bx + d g2 b 2 + d g3 b 3 ) 

+i 1 (d g0 b 1 - d gl b ) + i 2 {d g0 b 2 - d g2 b ) + i 3 {d g0 b 3 - d g3 b ) 
+ii(d g3 b 2 - d g2 b 3 ) + i 2 (d gl b 3 - d g3 bx) + i 3 (d g2 bi - d gl b 2 )} 
+ [-I a (d ea b + d e ibi + d e2 b 2 + d e3 b 3 ) 

+Ii(d e obi - deibo) + I 2 (deob 2 - d e2 b ) + I 3 (d e ob 3 - d e3 b ) 
+Ii(d e2 b 3 - d e3 b 2 ) + I 2 (d e3 bi - d e ib 3 ) + I 3 {d e \b 2 - d e2 b\)] 
+k b [Io(d g0 B a + d g iBi + d g2 B 2 + d g3 B 3 ) 

+Ii(d g0 B 1 - d gl B ) + I 2 {d g0 B 2 - d g2 B ) + I 3 {d g0 B 3 - d g3 B Q ) 
+Ii(d g2 B 3 - d g3 B 2 ) + J 2 (S s3j Bi - 8 gl B 3 ) + I 3 (d gl B 2 - d g2 B 1 )] 

+k b [{d e0 B + d el B x + d e2 B 2 + d e3 B 3 ) 

+i 1 (d el B Q - deo^i) + i2(d e2 B Q - d e0 B 2 ) + i 3 (d e3 B - d e0 B 3 ) 
+k{d e2 B 3 - 8 e3 B 2 ) + ^(d^i - el B 3 ) + i 3 (8 el B 2 - d e2 B x )\ 
-/i»(sg + hsl + i 2 s 9 2 + i 3 s 9 3 ) - /i«(I sg + hs\ + I 2 s e 2 + I 3 s e 3 ) 

-k bf i 9 (S 9 + hS? + i 2 S 9 + i 3 S§) - k b ^ e (I S c + hSt + I 2 S! + I 3 SI) (12) 

where, the first two terms arc the definition of momentum, and the last two terms are the definition of current; /u, 
/if, fi g and \i 9 are coefficients. 

The physical quantities in the quaternion space need to be extended to that in the octonion space. In the octonion 
space, the definition of physical quantity should be equal and consistent in the electromagnetic and gravitational 
interactions. 

The velocity Y(v , V\, v 2 , v 3 , V , V\, V 2 , V 3 ) of particle is defined as 

V = dR/dt = (v + i lV i + i 2 v 2 + i 3 v 3 ) + (I V + hVi + I 2 V 2 + I 3 V 3 ) (13) 

where, the first term is the velocity in G space, and the second term is the velocity in E space. 
In the octonion space, the charge and mass of particles are defined as follows 



So + *i s ! + *2 S 2 + *3S3 = Q 9 g (v + kvi + i 2 v 2 + i 3 v 3 ) (14a) 

Io4 + hst + I 2 s c 2 + I 3 s c 3 = Q g (I V + hV x + I 2 V 2 + I 3 V 3 ) (14b) 

S 9 + kS 9 + hSl + i 3 S§ = Q 9 (v a + iivi + i 2 v 2 + i 3 v 3 ) (14c) 

I S% + hSt + I 2 S e 2 + I 3 S e 3 = Qt(I V + hV! + I 2 V 2 + I 3 V 3 ) (14d) 



where, the first two equations are the definition of mass Q 9 g and Q g respectively; the last two equations are the 
definition of charge Qf and Q% respectively. The charges Q a g , Q g , Q 9 e and Q% can be known as the 'general charge' 
uniformly. 
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In the octonion space, there are four types of field sources in the electromagnetic-gravitational field. There exist 
two parts of field sources in the electromagnetic interaction, one part lies in G space (1, ^, ^3) and the other in E 
space (Jo, Ji, J2, la}- So does the gravitational interaction. 



C. Force and angular momentum 

In the electromagnetic-gravitational field, the force can be defined by the linear momentum P = pS/fi^, field 
strength B, and operator <>■ The field force F and variation rate J of field source are defined as 

F = E*oP, J=0*oP. (15) 

The definition of variation rate JJ of field source is different from that of field force F. But both of them have effects 
on the matter particle. The field force F includes the Lorentz force Q, while the variation rate JJ of field source 
includes the inertia force @. 

The angular momentum M is defined as (k rx is a coefficient) 

M = (R + k rx X) o P (16) 

The M is defined by the linear momentum P, field strength B, radius vector R, physical quantity X and operator 
0, and includes the angular momentum, work, kinetic energy, potential energy and force moment etc. 
The work and moment derived from the field force are defined uniformly as 

Wi = B o (R o P) + B o (k rx X o P) (17) 

where, the physical quantity B o (X o P) is a new part of the energy. 

The kinetic energy, potential energy and moment derived from the variation rate of field source are defined as 

W 2 = 0o(RoP) + 0o(fc„XoP) (18) 

where, the physical quantity (X o P) is the part of potential energy. 

Both definitions of energy Wi and W2 should exclude the potential energy part A o P, if the definition of angular 
momentum is confined in M = R o P . Hence this definition of angular momentum is inexact apparently. As shown in 
Aharonov-Bohm experiment [13] , the field potential A = o X has effects on the energy <0 o (X o P) , which includes 
the potential energy A o P . And the change of field potential can impact the movement of field source. 



IV. EQUATIONS SET OF ELECTROMAGNETIC-GRAVITATIONAL FIELD 

In the octonion electromagnetic-gravitational field, two types of the forces can be written together as follows from 
Eq.(15) (a — c is a coefficient) 

Z = F + aJ = (B + aO)*oP (19) 

And two types of energies can also be written as follows from Eqs.(17) and (18) 

W = Wi + aW 2 = (B + a0) o M (20) 

The unified definitions of force Eq.(19) and energy Eq.(20) show that, previous equations set of the quaternion 
electromagnetic field and octonion electromagnetic-gravitational field should be supplemented and generalized. It 
also shows that the octonion differential operator needs to be generalized to a new operator (0 + M/a). The 
physical characteristics of electromagnetic- gravitational field can be described uniformly from many aspects. 



A. Field source equation 

The extended definition of force Eq.(19) shows that, the field source § needs to be revised and generalized to the 
form of following equation from Eq.(12) 

fjS = -{M/a + 0)* o B = {fjSy - B* o M/a (21) 

where, (/xS)' = —0* B is the definition of the field source in Eq.(12) in the Maxwellian electromagnetic theory or 
the Newtonian gravitational theory. 

As one part of field source S, the term (B* o M/a) is directly proportional of field energy density (B* o M/ fJ&). The 
force-balance equation can be obtained when the force Z = from Eq.(19). And the variation of term (B* oM/a) will 
be one kind of necessary portion of force Z . 
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TABLE IV: The subficld types of electromagnetic-gravitational field. 



Operator 


Gravitational Interaction 


Electromagnetic Interaction 




pravitatioTi a l-pravitfltionfll sn hfiplH 


pi pet mm a pnpt.ic-pra vita tion a 1 sn HfiplH 


of G space 


G mass, Q 3 g 


G charge, Qf 




intermediate particle, 7^ 


intermediate particle, 7I 


operator Oe 
of E space 


gravitational-electromagnetic subfield 
E mass, Q g 

intermediate particle, jg 


electromagnetic-electromagnetic subfield 
E charge, Q% 
intermediate particle, 7! 



B. Power equation 

The extended definition of energy Eq.(20) shows that, the power N needs to be revised and generalized to the form 
of the following equation. 

N=(B + aO)*oW (22) 

In the electromagnetic-gravitational field, the physical quantity X has effect on field potential X , angular 
momentum X o P , energy Bo(XoP) and power (BoB')o(XoP) etc. The introduction of physical quantity X makes 
the definition of angular momentum M and energy W more integrated, and the theory more self-consistent. 

In Eq.(20), the conservation of angular momentum in electromagnetic-gravitational field can be gained when W = 0. 
And the energy conservation of electromagnetic-gravitational field can be attained from Eq.(22) when N = 0. 

C. Dark matter 

In the octonion spacetime, there exist four types of subfields and their field source of electromagnetic-gravitational 
field. In the Table 4, the electromagnetic-gravitational (E-G) subfield and gravitational-gravitational (G-G) subficld 
are 'electromagnetic field' and 'modified gravitational filed' respectively. Their general charges are G charge and 
G mass respectively. The electromagnetic-electromagnetic (E-E) subfield and gravitational-electromagnetic (G-E) 
subfield are both long range fields and candidates of the 'dark matter field'. Their general charges (E charge and E 
mass) are candidates of 'dark matter'. The physical features of the dark matter meet the requirement of Eqs.(10), 
(11), (16), (19), (20), (21) and (22) . 

In the general charge (Q 9 , Q g , Q 9 and Q%) and intermediate particle (t|, 7?, 7| and two types of general 
charges (familiar charge Q 9 and mass Q 9 g ) and one type of intermediate particle (photon j 9 ) have been found. Hence 
the remaining two types of general charges and three types of intermediate particles are left to be found in Table 4. 
From Eq.(12), we find that the electromagnetic-gravitational field can transform one form of intermediate particle 
into another. In other words, the photon j 9 can be transformed into the j 9 , 7|, or 7! . 

The particles of ordinary matter (electron and proton etc.) possess the G charge together with G mass. The 
particles of dark matter may possess the E charge with E mass, or G mass with E charge, etc. It can be predicted 
that field strength of electromagnetic-electromagnetic and gravitational-electromagnetic subfields must be very weak 
and a little less than that of the gravitational-gravitational subficld. Otherwise they should been detected for a long 
time. Hence the field strength of electromagnetic-electromagnetic and gravitational-electromagnetic subfields may be 
equal, and a little less than that of gravitational-gravitational subfield. 



TABLE V: The comparison between ordinary matter with dark matter. 





Ordinary Matter 


Dark Matter 




subfield 


G-G subfield 


E-G subfield 


G-E subfield 


E-E subfield 


field potential 


OgOXg 


Og O X e 


Oe OXg 


Oe O X e 


field strength 


OgOAg 


9 A e 


Oe O A 9 


Oe O A e 


field source 


9 °B 9 


0; oB e 


0: oI 9 


Oe °Be 
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FIG. 2: Linear movement. The change of spatial distribution of mass m will affect free-fall movement of particle m. 
The linear movement status of particle m in point a and b are different, though the field strength are equivalent in 
those points. 



V. EXPERIMENTS AND PHENOMENA 



The field sources (current and momentum) in Eq.(19) stays in the force-balance status in the electromagnetic- 
gravitational field when Z = 0. The influence between the 'dark matter field' (electromagnetic-electromagnetic 
and gravitational-electromagnetic subfields) and the 'modified gravitational field' (gravitational-gravitational sub- 
field) on the force Z can be discussed with this feature. The conservation of angular momentum in Eq.(20) in the 
electromagnetic-gravitational field can be achieved when W = 0. The influence between the 'modified gravitational 
field' and 'dark matter field' on the angular momentum and energy can be considered by this property. It is simpler 
and more convenient to represent the field source (fjE)' as /jS> in the following discussion. 



A. Linear movement 



In the gravitational-gravitational subfield (<?oi ; 0, 0) in Fig. 2, the particle m is in free-fall movement and its momen- 
tum is (sg, sf , 0, 0). When the particle m is in weightlessness status, 

Z = (nsoi + aO)* o { M |( s g + i lS {) - glJa) / M » = (23) 

then eight equations can be obtained from Eq.(23). And its equation in i\ direction is 

Msfoisg + otfj, 9 g (d gl s% - d g0 s g ) - (<9 9 i3oi + 5oi/«) = ( 24 ) 

When the last term of the above equation is zero, we can achieve 

9ois 9 Q = c(d g0 s 9 - d g isl) (25) 

where, the physical quantity can be written in familiar form, c/oi — #/c , Sq = mc , s 9 = mv\. 

The above shows that, when the spatial distribution of the mass m changes, as one part of inertia force, <9 9 iSq will 
affect free-fall movement of particle m. The omitted term (<9 fl igoi +5oi/ a ) m Eq.(24) will influence free-fall movement 
of particle m also. Then the linear movement status of particle m in the point a and b are different in Fig. 2, though 
the field strength goi in the point a and b are equal. 



B. Stationary status 



In the area, where coexist the electromagnetic-gravitational subfield (fooij0,0) and gravitational-gravitational sub- 
field (goijO,0). The electric charged particle q is in the stationary status in the area in Fig. 3, its momentum is 
(sq, 0, 0, 0) and its current is (Sq , 0, 0, 0). When the particle q stays in the force-balance status 

Z = {(k9oi + fcaii&oi) + «0}* o {(n g s a Q + k b ^I S 9 ) - (g 2 01 + k 2 a b 2 01 )/a} - (26) 

so its equation in i\ direction is 

(^9014 ~ k bf iik a b Q1 S 9 ) + a(^ g d gl s 3 Q - k bf ild el S 9 ) - d gl (g 2 Q1 + k 2 a b 2 Q1 ) - g 01 (gl + k 2 a b 2 Q1 )/a = (27) 
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FIG. 3: Stationary status. The transform of energy density in the electromagnetic-gravitational field will affect the 
force-balance of charged particle including the stationary status of electric charged particle. In consequence the 
stationary status of electric charged particle q in point a and b will be different. 



The last two terms of the above equation show that, the transform of energy density in the electromagnetic- 
gravitational and the gravitational-gravitational subfields will affect the force-balance of charged particle. When the 
sum of above last three terms equals to zero, and k a k b — [i 9 g j [i 9 e , we can conclude from above 

5oiSo ~ koi^o = 

where, the physical quantity can be written in the familiar form, 6qi = E/c, 501 = g/c, Sq — mc, Sq = qc . 

The last three terms in Eq.(27) will impact the stationary status of electric charged particle. So the stationary 
status of electric charged particle q in point a and b will be different in Fig. 3, though the field strength 601 in point 
a and b are equal. 



C. Planar circumference movement of dark matter 



In the electromagnetic-gravitational subfield (0, 0, &12), the charged particle q is in the planar circumference 
movement in Fig. 4. Its momentum is (sq, sf , s 9 , 0), while its two types of currents are (S§, Sf,S%, 0) and (Sq, fif , Sf, 0) 
respectively. When the particle q is in the force-balance status, 

Z = (I 3 k a b 12 + aO)* o 

{^(s 9 + iisf + hs 9 2 ) + hvlihSl + hS{ + I 2 S C 2 ) + k bf i 9 (S 9 + nSf + hSi) - k 2 a b 2 12 /a} j pP g = (28) 
so its equations in i\ and i 2 directions are respectively 

kahbunlSZ + an B g (-d g os{ + <9 g isg - d g3 s^) - d g i(k a bi 2 ) 2 

+ak bf i e e (d e0 Sl - d el S e G + d e3 S e 2 ) + ak b ^~d g oSf - d gl S$ + d g3 S 9 2 ) = (29a) 
~k a k b b 12 n L e Si + a^-dgosl + d g2 s 9 + d g3 s{) - d g2 (k a b 12 ) 2 

+ak bf i e e (d e0 S e 2 - d e2 S e G - d e3 SD + ak b ^ e (-d g0 S 9 + g2 S 9 + 8 g3 S 9 ) = (29b) 

It makes out that field source (Sq , Sf, S%, 0) of 'dark matter field' has direct influence on the movement of field 
sources (Sq, Sf, S 2 ,0) of familiar 'electromagnetic field'. 

When the sum of last three terms in above each equation equals to 0, and k a k b = H 9 g / H 9 e , the force-balance equations 
of charged particle, which are in the circumference movement in electromagnetic field and modified gravitational field, 
can be achieved from above 

(J4HM )b 12 S! + c(-d g0 s 9 + d gl s 9 - 8 g3 s 9 2 ) = (30a) 
(n 9 g H e Jl4 )b 12 Sf - c(-d gQ 4 + dg24 + d g3 s 9 ) = (30b) 

The above equations show that, as the parts of inertia force, the terms (dgiSy — d g3 s 2 ) and (d g2 s 9 Q + d g3 sf) will 
affect the planar circumference movement of charged particle. The terms d g ib 2 2 in Eq.(29a) and d g2 b\ 2 in Eq.(29b) 
will affect the planar circumference movement of charged particle q also. 
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FIG. 4: Planar circumference movement. The 
field source of 'dark matter field' has direct 
influence on the movement of field source of 
familiar 'electromagnetic field'. 



FIG. 5: Near circularity and rotation. The omitted 
'magnetic' part of gravitational-gravitational subfield 
strength affects the movement of particle, and is one of 
important factors resulting in rotation of celestial bodies. 



D. Near circularity and rotation 

The large-mass object M generates the field strength (301,302,512) of the gravitational-gravitational subfield in its 
equator plane. The particle m moves around M in planar helix, and its momentum in orbit plane is (sq, sf , s 2 , 0). 
The orbit plane of particle m is superposed with the equator plane of large-mass object M in Fig. 5. When particle 
m stays in force-balance status, 

Z = (*i5oi + ^2.902 + «3.9i2 + a0)* ° {a^Oo + *i s i + ^ s l) ~ B V a } I V g g 
so its equations in i\ and i 2 directions are 

A*g(ffoiSo + 912S2) + a V 9 g (-dgO s i + d gi s o ~ 9g34) - d gi B2 - 9oiB 2 /a 
Vg(go2S 9 - gi2s{) + afi 9 g (-d g0 s?, + d g2 s 9 + d g? ,s{ ) - d g2 B 2 - g 02 B 2 /a 

where, B 2 = g 2 x + g 2 2 + g\ 2 . 

When the sum of last two terms in above each equation is equal approximately to zero, we can attain 

(301S0 + S12S2) + ci-dgosf + dgis 9 - d g3 s 9 2 ) = (33a) 
(902S0 - 012 *f) + c(-d g0 s 9 2 + d g2 s 9 + dg 3 sf) = (33b) 

When 312 = 0, we can learn that the particle m moves around object M in planar ellipse from above equations. 
That is the near circularity of planets. When 301 = 302 = 0, it can be found that the particle m moves in planar 
circumference from above equations. That is the rotation of planets. Therefore, in gravitational-gravitational subfield 
(301,302,312) , the centroid movement of particle m is the superposition of these two types of movements. That is, 
the center of planar circumference movement of particle moves in planar ellipse. When the particle m moves around 
big mass M, the movement of particle m is in revolution, rotation and swing because of weakness and nonuniform 
distribution of 312 [11]. It illustrates that, the omitted 'magnetic' part 312 of gravitational-gravitational subfield 
strength affects the movement of planet m, and is one of important factors resulting in rotation of celestial bodies. 

As shown above, the terms (<9 g iSg — d g 3S 2 ) and (dg 2 SQ + 9 g 3sf), as certain parts of inertia force, will affect the 
planar helix movement of particle in Fig. 6. In the solar system, when spatial distribution of planetary mass and 
field energy density are invariable, or when the planetary momentum is invariable along the direction perpendicular 
to the Ecliptic Plane, the movement of planet is consistent with the Newtonian gravitational theory. Contrarily, 
in the Milky Way Galaxy, when the spatial distribution of solar system mass and field energy density are variable, 
or when the angular momentum of the solar system is variable along the direction perpendicular to the Galactic 
Equator Plane, the centrifugal force of solar system will change, and then the movement of solar system will betray 
the familiar gravitational theory. If the change of mass, energy density and momentum lead to counteract the certain 
parts of centrifugal force, it will result in the revolution speed increase to generate enough centrifugal force to balance 
gravitational interaction of the Milky Way Galaxy. Therefore the inferences may solve inconsistency between the 
overspeed rotation of galaxy and deficiency of gravitation in some extent [13]. 



= 

= 
= 



(31) 

(32a) 
(32b) 



E. Near coplanarity and corevolving 



On the orbit plane of particle m which momentum is (sg, s 9 , sf , sf), the gravitational-gravitational subfield (301 + 
323, 302 + 331, 303 + 312) and the electromagnetic-gravitational subfield (6qi + ^23, &02 + &3i, &03 + &12) are generated by 
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FIG. 6: Overspeed revolution. If the change of mass, 
energy density and momentum lead to counteract the 
certain parts of centrifugal force, it will result in the 
revolution speed increase to generate enough 
centrifugal force to balance gravitational interaction. 




FIG. 7: Near coplanarity and corevolving. When the 
sense of revolution of particle m is in the same sense 
of rotation of object M, the movement of particle m is 
attracted by the force which points to the orbit plane, 
hence the movement of particle m is stable. 



the large-mass object M which owns rotation and charge in Fig. 7. When the particle m, which moves around M in 
spatial helix, is in the force-balance status, 

Z = [h(goi + .923) + h(902 + 931 ) + *3(.903 + 512) 

+k a {Ji(&oi + b 23 ) + h(b 2 + 631) + h(b 3 + 612)} + «<>]* 

{/**(*§ + i lS { + hsl + i 3 s 9 ) - B 2 /a) ln% = (34) 
so its equations in i\ , i 2 and ^3 directions respectively are, 

H 9 g {(m + 923)4 + (.902 + .931)4 - (503 + 312)4} 
+an 9 g (-d g Qs{ + d gl s 9 + d g2 4 ~ d g3 s 9 2 ) ~ d g iB 2 - (g 01 + g 23 )B 2 /a = 
^g {(302 + 33i)4 ~ (.9oi + .923)4 + (303 + 312)41 

+a f i 9 (-d g0 s 9 - d g i4 + d 9 24 + d g3 4) ~ ^ 2 S 2 - (. 90 2 + 9 3 i)B 2 /a = 

{(.903 + .912 ) 4 + (301 + 323 ) 4 - (302 + 33l)4} 

+aii g (-d g0 s 9 3 + d gl s 9 - d g2 s 9 + d g3 s 9 ) - d g3 B 2 - (g 03 + 9l 2)B 2 /a = 

where, B 2 = (g 01 + g 23 ) 2 + (.902 + 33i) 2 + (303 + 312) 2 + k 2 a {(b 01 + b 23 ) 2 + (602 + 631) 2 + (b 3 + h 2 ) 2 }- 

The last two terms of above equations show that, the change of energy density (S 2 //x|) in the gravitational- 
gravitational and electromagnetic-gravitational subfields will influence the movement of particle m. When the sum of 
last two terms in above each equation is zero, the equations of gravitational-gravitational subfield can be summarized 
as follows 

{(301 + 323 ) 4 + (.902 + 33i)4 _ (.903 + 312)4} + c (-9go4 + d gi4 + d g^4 ~ ®gz4) = ( 36a ) 
{(302 + 33i)4 ~ (.9oi + 323)4 + (.903 + 3i2)s?} + c(-d g o4 + d gz4 + d gs4 ~ d gi4) = ( 36b ) 
{(303 + 312)4 + (.901 + 323)4 - (302 + 33l)4} + c (-9go4 + d a^4 + d g^4 ~ d g?4) = ( 36c ) 

From above equations we can learn that, when the sense of revolution of particle m is in the same sense of rotation 
of object M, the movement of particle m is attracted by the force which points to the orbit plane, hence the movement 
of particle m is stable. That is the near coplanarity of planets. The particle m, which moves around object M in 
spatial helix, revolves and waves up and down slowly around the equator plane of object M (the earth's movement 
relatives to the equator plane of the sun; the movement of sola r sy stem and pulsar relative to the galactic plane of 
the Milky Way Galaxy) , at the same time it rotates with swing hj. When the sense of revolution of particle m is in 
opposite sense of rotation of object M, the movement of particle m is repelled by the force which deviates from the 
orbit plane and hence the movement of particle m is unstable. It shows that, the sense of revolution of particle m 
affects directly its movement stability. That is corevolving (or prograde) of planets. Therefore, the solar system and 
galaxies are inclined to thin disk structure as a whole [3] • 

To be extended, two galaxies rotating in the same sense of rotation are in attraction and revolution around each 
other, even shorten their interval and swallow up each other to become one galaxy. Two galaxies rotating in opposite 



(35a) 
(35b) 
(35c) 
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FIG. 8: Spatial helix movement. In Eq.(39), when the 
right side of equal mark in last equation is equal to 
zero, and the first three equations are not all zero, the 
angular momentum M| is close to be in conservation 
and in consequence more complicated situations will 
appear on the revolution orbit of particle m. 



FIG. 9: Planar helix movement of celestial body. It 
shows that there are the coupling influences between 
the 'dark matter field' and 'modified gravitational 
field'. When the last term is zero in Eq.(42), there 
will be no influence of the 'dark matter field' upon the 
'modified gravitational field'. 



senses of rotation will repel and increase their interval. Two galaxies rotating in arbitrary senses of rotation will be 
affected by more complicated forces and moment [TEj . 



F. Spatial helix movement 



On the orbit plane of particle m, which angular momentum is (Mq , Mf , M|, M3), the gravitational-gravitational 
subfield (501 + S23 > 902 + <?3i > <7o3 + 912) is generated by large-mass object M which maintain its rotation simultaneously 
in Fig. 8. When the particle m, which moves around M in spatial helix, stays in conservation of angular momentum, 



W = {ii(sox + 323) + h{go 2 + 331) + ^3(503 + .912) + aO} ° (M$ + i\M{ + i 2 Mf 







so its equations in 1, and i 3 directions are respectively 



c(d g0 M% - d gl M? - d g2 Ml - d g3 M§) - {(901 
c(d g0 Ml + d gl Ml + d g2 M§ - d g3 M!) + {(301 
c(d g0 M$ - d B iM§ + d g2 M$ + d g3 M°) + {(302 
c(d g0 M§ + d gl Ml - d g2 M° + d g3 M$) + {(303 



92 3 )M( + (302 + g 3 i)M$ + (303 
92 3 )Mg + (302 + g 3 % )Mf - (303 
33i)M(? - (301 + 92 3 )M§ + (303 
S12 )M,f + (301 + 92 3 )Ml - (302 



9i2)M§} 
9i2)M%} 

<?12)Mf} 

33i)Mf} 



(37) 

(38a) 
(38b) 
(38c) 
(38d) 



In the gravitational-gravitational subfield, above equations set show that, the conservation of angular momentum 
will be affected by the field strength, angular momentum and its variation rate etc. When Mf — A j9 — 
obtain from the above equations set 



M% = 0, we can 



d, 



d g2 M§ 
d gl M§ 
d g0 M§ 



d g0 M°-{(g 03 + 
-d gl M$ - {( 501 
d„2M° + {{g 02 + g 31 )Ml 



-d, 



{(50 



9i2)M§}/c 
+" 92 3 )Ml - (302 

- (501 + 
5i2)M 9 }/c 



+ 33i)M 3 9 l/c 
92s)M§}/c 



(39a) 
(39b) 
(39c) 
(39d) 



When all terms in the right side of equal marks of the above equations set are zero, the angular momentum A/| of 
particle m along the sense of revolution is in conservation. And if the right side of equal mark of the last equation 
is equal to zero approximately, whilst the right side of equal marks of the first three equations are not all zero, the 
angular momentum M| is close to be in conservation and consequentially more complicated situations will appear on 
the revolution orbit of particle m. For example, the asteroid will probably come near to the planet gradually, and the 
moon may perhaps keep away from the earth in the solar system [ltH Il7| . 



G. Planar helix movement of dark matter 



The field strength of gravitational-gravitational and gravitational-electromagnetic subfields in Fig. 9, which in equa- 
tor plane of large-mass object M, are (301,302,312) and (Gqi, G02, G\ 2 ) respectively. The particle m moves around 



13 



M in planar helix, and its two types of momentums in orbit plane are (sq, sf, s 2 ,0) and (sq, s\, s 2 ,0) respectively. 
The orbit plane of particle m is superposed with equator plane of large-mass object M. When the particle m is in 
force-balance status, 

Z = (iiffoi + *2.9o2 + kg\2 + iiGoi + hG 02 + J3G12 + aO)* o 

{/i»(ag + iisf + i 2 sf ) + n' g {hs% + list + hsl) - B 2 /a} /n° = (40) 

so its equations in ii and i 2 directions are respectively 

Mg(foiSo - 9i2s 9 2 ) + a^ 9 g {-d g os{ + d g i4 - d g3 s 9 2 ) 

+ M g {a(d e o4 - deiso + d e3 s e 2 ) + (G 12 s e 2 - G is e )} - {g 01 B 2 /a + d gl B 2 ) = (41a) 

+H e g {a(d e0 4 - d e2 s e - d e3 st) - (G 12 s\ + G 02 s e )} - (g 2B 2 /a + d g2 B 2 ) = (41b) 

where, B 2 = g 2 01 + g 2 2 + g 2 2 + G 2 01 + G 2 2 + G\ 2 . 

When the last term of above each equation is zero approximately, we can conclude 

(5oi«o - S12S2) + ci-dgosf + d gl s 9 - d g3 s 9 2 ) 

+ 04/^) {c(d e0 st - 8 el s e + d e3 s e 2 ) + (G 12 s e 2 - Goi^)} = (42a) 
(502S0 + ffi2sf) + c{-d g0 s?, + d g2 s 9 Q + 9 g3 s?) 

Hl^V^) {c(d e0 s e 2 - d e2 sl - d e3 sl) - (Giasf + G 02 s e )} = (42b) 

In above equations set, the above equations are the movement equations of gravitational-gravitational subficld, 
including the affecting from the gravitational-electromagnetic subfield. It figures out that there are some coupling 
influences between the 'dark matter field' (gravitational-electromagnetic subfield) and 'modified gravitational field' 
(gravitational-gravitational subfield) . When the last term is zero in the above each equation, there will be no influence 
of the 'dark matter field' (gravitational-electromagnetic subfield) upon the 'modified gravitational field' (gravitational- 
gravitational subfield). 

If the preceding field strength and field source are extended to four types of field strength and field source of 
electromagnetic-gravitational field respectively, many extra and also more complicated equations set of interplay 
between 'dark matter field' and 'modified gravitational field', will be gained from Eqs.(19) and (20). 

VI. CONCLUSIONS 

In the observed phenomena of the celestial body, which are inconsistent with the results of current gravitational 
theory, some are caused by the modified gravitational field (gravitational-gravitational subfield) in the electromagnetic- 
gravitational field, and some are caused by the dark matter field (electromagnetic-electromagnetic and gravitational- 
electromagnetic subficlds) . 

The gravitational-gravitational subfield is the modified gravitational field, which includes familiar Newtonian gravi- 
tational field. The features of modified gravitational field predict that, (1) the planetary orbits possess near coplanarity, 
near circularity and corevolving, and the planets own rotation property, (2) the centrifugal force of celestial body will 
change and lead to fluctuation of revolution speed, when the field energy density or the angular momentum of celestial 
body along the sense of revolution varied, (3) the two galaxies rotating in the same sense of rotation will attract and 
revolve around each other, hence shorten their interval and swallow up each other to become one galaxy. 

The electromagnetic-electromagnetic and gravitational-electromagnetic sub-fields are both long range fields and 
candidates of dark matter field. Their general charges are candidates of dark matter. The combination of above general 
charges and the general charges of gravitational-gravitational or electromagnetic-gravitational subfields are candidates 
of dark matter also. The field strength of electromagnetic-electromagnetic and gravitational-electromagnetic subficlds 
may be equal and very weak, and a little less than that of gravitational-gravitational subficld. The features of dark 
matter predict that, (1) the dark matter field possesses two sorts of general charges makes the dark matter particles 
very diversiform, (2) the field strength of two types of dark matter fields are both weaker than that of the Newtonian 
gravitational filed, (3) the distribution of dark matter will affect the stabilizing velocity of celestial body. If the 
field sources of the electromagnetic-electromagnetic and gravitational-electromagnetic subfields distribute properly, 
the obtained gravity of celestial body will increase, and the stabilizing velocity of celestial body will also increase 
correspondingly. 
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